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Progress using generalized lattice Dirac operators to parametrize the 
Fixed-Point QCD action 
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We report on an ongoing project to parametrize the Fixed-Point Dirac operator for massless quarks, using 
a very general construction which has arbitrarily many fermion offsets and gauge paths, the complete Clifford 
algebra and satisfies all required symmetries. Optimizing a specific construction with hypercubic fermion offsets, 
we present some preliminary results. 



1. Motivation for improved fermionic ac- 
tions 

There are many reasons to seek an improved 
lattice Dirac operator. The most commonly used 
one, the Wilson Dirac operator, has a number 
of problems. Chiral symmetry, an essential el- 
ement of QCD, is broken explicitly by a term 
which removes the fermion-doublers and so the 
pion is not massless for zero quark mass. The 
pion can be tuned massless by varying the hop- 
ping parameter. However, this does not remove 
other problems, such as the mixing of operators 
in different chiral representations, which makes 
the extraction of chiral matrix elements difficult. 
There are also problems measuring chiral quanti- 
ties such as the quark condensate and zero modes 
of the Dirac operator and connecting these to the 
topology of the gauge fields, due to, among other 
things, topological defects and exceptional config- 
urations. In addition, the Wilson Dirac operator 
has large lattice artifacts. 

A solution to problems of chirality has been 
proposed by Ginsparg and Wilson H: a lattice 
Dirac operator which satisfies 



{ 75 ,C- 1 } = 2 fl i? 75 , 
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where R is any local function which commutes 
with 75 and a is the lattice spacing, was expected 
to define a lattice theory with the same chiral 

*Based on a tafk by K.H. and a poster by T.J. 



properties as the continuum theory. Unfortu- 
nately, no solution to this relation was found and 
the idea was abandoned. The first lattice regu- 
larization of fermions with chiral symmetry, the 
Domain Wall fermions Q and the related Over- 
lap construction |3| , followed a different path and 
seemed to be unrelated to the Ginsparg- Wilson 
(GW) relation. 

Following the observation that the Fixed-Point 
Dirac operator satisfies the GW relation Q, the 
interest turned to this general formulation again. 
The GW relation turned out to be a powerful 
theoretical tool: it implies an index theorem on 
the lattice || and all the chiral properties of 
the formal continuum theory ||. This last fea- 
ture became trivial after it was recognized that 
an exact chiral symmetry transformation follows 
from cq.(|l|) 0. The observation that the Over- 
lap Dirac operator satisfies the GW relation con- 
nected the Domain Wall approach to this gen- 
eral formulation |J| . The progress in vector chiral 
symmetry was followed by a breakthrough in the 
formulation of regularized chiral gauge theories 
also §. 

All these developments have led to much recent 
activity on the theoretical understanding of lat- 
tice chiral symmetry pp| , pl| and on constructing 
improved lattice Dirac operators which approxi- 
mately satisfy the GW relation At the 
same time, detailed tests and the first large-scale 
simulations with the simplest version of Domain 
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Wall and Overlap Dirac operators demonstrated 
the power and difficulties of working with chirally 
symmetric QCD actions [lq| . 

2. General lattice Dirac operators 

Before we explain our approach to the 
parametrization of the Fixed-Point Dirac oper- 
ator we briefly discuss the structure of general 
Dirac operators which satisfy the basic symme- 
tries on the lattice. This will make clear what 
kind of lattice operators can be used in the 
parametrization of the Fixed-Point Dirac oper- 
ator or in any other parametrization of a Dirac 
operator. A more detailed discussion of these is- 
sues is given in pifl . A general gauge covariant 
lattice operator with color, space and Dirac in- 
dices can be written as: 

D=J2T A J2c(T A ,l)U(l), (2) 

A I 

where the Ta are elements of the Clifford algebra 
basis which we define as T = 1, 7^, ?<7 M „, 75, 7^75, 
I denotes a path, U(l) the product of the link ma- 
trices along this path and c(Ta, I) is the coupling 
for the given path and Clifford algebra element. 

The basic symmetries of the Dirac operator, 
which are gauge symmetry, 75-hermiticity, charge 
conjugation, hypercubic rotations and reflections, 
impose the following restrictions on eq.(^J). 

• Translation invariance requires that the 
Dirac operator D(n, n + r) connecting lat- 
tice sites n and n + r depends on n only 
through the n-dependence of the gauge 
field, i.e. the couplings c(Ta,1) are con- 
stants or gauge invariant functions of gauge 
fields, respecting locality and invariance un- 
der the symmetry transformations. 

• 75-hermiticity and charge conjugation to- 
gether imply that the couplings c(Ta, I) for 
our choice of the Clifford algebra basis are 
real. Further, from hermiticity it follows 
that the path I and the opposite path I (or 
equivalently, U(l) and 17(1)^) should enter 
in the combination 

r(u(l) + e r U(iy) , (3) 



where the sign er is defined by 7 5 r^7 5 = 
erT. 

• Permutations and reflections of the coordi- 
nate axes (hypercubic rotations) imply that 
for a given path Iq a whole class of paths 
belongs to the Dirac operator. These paths 
are related to Iq by all the 16 x 24 = 384 
reflections and permutations of the coordi- 
nate axes. Under such a symmetry trans- 
formation a — 1, . . . , 384 the Clifford alge- 
bra element Tq associated with Iq generally 
transforms into r' tt ). Furthermore the sign 
of the couplings may change, whereas their 
absolute value remains unchanged. 

A Dirac operator satisfying all the basic sym- 
metries can be written as 

d = ]r C (r ,ME° Q 

O a = T<- a) (d(l^) + e r £>(2 ( " ) ) t ) ■ (4) 

where the sum runs over a set of reference paths 
defined by Tq and Iq as well as over all the symme- 
try transformations a defined by the group of per- 
mutations and reflections of the coordinate axes. 

It is quite natural to include the full Clifford 
algebra in the Dirac operator. The Wilson opera- 
tor already contains 1 and 7^, the 0(a) Symanzik 
condition involves a^ v and the topological charge 
is proportional to Tr^D), which is zero unless D 
contains 75. Furthermore, for 75-hermitian Dirac 
operators, the Ginsparg- Wilson relation ([j]) can 
be written as 

D + = 2 a D^RD, (5) 

which can only be satisfied if D contains the com- 
plete Clifford algebra. 

3. Efficient implementation of general 
Dirac operators 

On first appearances one might think that it is 
not feasible to calculate such a general Dirac op- 
erator which has many different couplings, where 
every coupling might contain as many as 768 
paths. But one has to keep in mind that the cal- 
culation of propagators for small quark masses 
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needs several hundreds or even thousands of con- 
jugate gradient steps and therefore one can afford 
to spend some time to precalculate and store the 
whole Dirac operator before one starts the calcu- 
lation of the propagator. On top of this there are 
two reasons why the building of general operators 
can be done in a very efficient way: 

• There are a lot of paths which are invariant 
under certain subgroups of the reflections 
and permutations which reduces the num- 
ber of terms significantly and in some case 
even leads to a cancellation of certain terms 
because they have relative minus signs. 

• A less trivial fact is that the sum of paths 
for many couplings can be factorized in an 
efficient way, i.e. that such large sums of 
many paths can be written as a product of 
smaller sums of fewer paths. 

As an example we consider a nearest neighbour 
coupling with Tq = 75 and generating path Iq = 
[2, 1, —2, 3, 4, —3, —4], where a path I is a sequence 
of steps in the ±/z direction. All the paths of this 
coupling can be written in the following compact 
way: 

75 ^ tp,vpa {SftuPpa + Pp a Sp, v + h.C.) , (6) 

where the color matrices and P pa are cer- 
tain combinations of staples or plaquettes, respec- 
tively. When all the plaquettes and staples and 
the most frequent combinations, e.g. Ppa, are pre- 
calculated then operators as in eq.(^) can be cal- 
culated very quickly. As an illustration of this 
we consider one of our parametrizations of the 
Fixed-Point Dirac operator. It has 39 couplings 
in total, at least one per offset on the 1-hypercube 
and per type of Clifford algebra element. Build- 
ing this operator takes roughly 15 times as long as 
to multiply the operator with a vector and there- 
fore it is a very small fraction of the time used to 
perform a calculation of one propagator. 

The basic operation required to calculate prop- 
agators or eigenvalues is the multiplication of the 
Dirac operator times a vector. For a Dirac oper- 
ator with hypercubic fermion offsets which con- 
tains the complete Clifford algebra, the matrix 



times vector multiplication is roughly 40 times 
more expensive than for the Wilson Dirac opera- 
tor. 

In p4| , we have included a MAPLE code which, 
given an input path I and any Clifford algebra 
element T, gives all of the paths l^ a ' and algebra 
elements generated by all permutations and 
reflections of the coordinate axes. This can be 
used by anyone who wishes to use a very general 
lattice Dirac operator, independently of how this 
operator is optimized. 

4. Fixed-Point improved actions 

The Fixed-Point method of improving lattice 
actions is inspired by the Rcnormalization Group 
flow of asymptotically-free theories Jll| . Consider 
a lattice action which contains all possible inter- 
actions. The Renormalization Group is defined 
by some blocking function which averages over 
the fields to produce a new action on a coarser lat- 
tice with fewer fields. The new action typically 
has different couplings from the original action 
and so we can imagine the blocking step as a flow 
in the coupling parameter space. Infinitely many 
blocking steps give the complete Renormalization 
Group (RG) trajectory. In Figure [|, we show 
the RG trajectory for QCD with massless quarks. 
The Fixed-Point (FP) has the property that the 
couplings are reproduced after a blocking step. 
For asymptotically-free theories, the Fixed-Point 
is on the surface of vanishing coupling g = 0. 

If one starts on this surface, the RG trajectory 
flows quickly to the Fixed-Point. If one starts 
close to this surface at some small coupling g 
(i.e. small lattice spacing a), the RG trajectory 
flows quickly towards the Fixed-Point and then 
flows away from it. Starting on the RG trajec- 
tory with arbitrarily fine lattices with arbitrarily 
small lattice artifacts, one can reach any point on 
the trajectory by making sufficiently many block- 
ing steps. All actions on the trajectory describe 
the same physics. The physical observables of the 
continuum quantum theory are identical to those 
of any lattice quantum theory on the RG tra- 
jectory, independently of the lattice coarseness. 
Such lattice actions are said to be quantum per- 
fect. The Fixed-Point action is an approximation 
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Figure 1. Rcnormalization Group trajectory of 
asymptotically- free theories. 



to the RG trajectory for small couplings g and is 
classically perfect, i.e. it completely describes the 
continuum classical theory without discretization 
errors. 

Fixed-Point actions have many desirable fea- 
tures. By closely approximating the RG trajec- 
tory, they are expected to have much reduced 
quantum lattice artifacts. They can be optimized 
for locality. The Fixed-Point Dirac operator satis- 
fies the Ginsparg- Wilson relation and so has good 
chiral behavior. The Fixed-Point QCD action has 
well-defined topology and satisfies the index the- 
orem on the lattice. Although more expensive 
than standard actions, studies of other models 
have shown that practical Fixed-Point actions can 
be constructed and used [0. 

5. Free Fixed-Point fermions 

We first consider free massless fermions. Be- 
cause the fermionic action is quadratic in the 
fermionic fields, the Renormalization Group step 
for the fermionic fields amounts to Gaussian inte- 
gration, which can be done exactly. The blocking 
step connects the Dirac operator on the fine (F) 
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Figure 2. Free fermion dispersion relation. 



and coarse (C) lattices by 

Dc 1 = - + D F 1 n\ (7) 

provided Dp has no zero modes, where k is an 
optimizable free parameter of the blocking and £1 
is the blocking function used to integrate out the 
fine fields to produce the coarse fields. The FP 
Dirac operator is reproduced under the blocking 
step i.e. Dc = Dp = Dpp. For free fermions, this 
equation can be solved exactly analytically. The 
FP Dirac operator is local, i.e. the couplings fall 
off exponentially with distance, and the rate of 
fall off can be maximized by varying n. However, 
Dpp contains infinitely many couplings. For prac- 
ticality, Dfp is approximated with an ultra-local 
operator, for which each point is only coupled to 
its neighbours on the hypercube. In Figure ^, we 
see the exact Dfp reproduces exactly the contin- 
uum dispersion relation, while the approximate 
FP Dirac operator has much smaller discretiza- 
tion errors than the Wilson operator. The un- 
physical branches in the approximate -Dfp occur 
at large momenta and only have a small effect. 
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Figure 3. Free fermion eigenvalue spectrum of the 
truncated FP Dirac operator. 



The Ginsparg- Wilson relation eq.(0) is a con- 
straint on {75, D -1 }. The FP equation (Q) con- 
nects the propagators D^ 1 on the coarse and fine 
lattices. Combining these equations gives 



r c = - + n R F fl\ 

K 



(8) 



and at the Fixed-Point, Rc = Rf = Rfp- For 
free fermions, this equation can also be solved 
exactly analytically. Choosing a symmetric over- 
lapping block transformation Q with a scale factor 
2 averaging over hypercubes |lq] , the exact Rfp 
is ultra-local and has only hypercubic couplings 
— no approximation for practicality is required. 
The block transformation determines R. With 
other methods to build a Dirac operator satis- 
fying the Ginsparg- Wilson relation, for example 
the Overlap construction, R is unconstrained and 
typically R = |. 

Defining a rescaled Dirac operator d — 
\2R D \f2R, the Ginsparg- Wilson relation can 
be written (with lattice spacing a = 1) 



i.e. the eigenvalues of d lie on a circle of radius 
1 centred at (1,0). Using Rpp and the hypercu- 
bic approximation of Dpp , we show the eigenvalue 
spectrum of d in Figure ||. It lies almostly exactly 
on the circle, indicating that the hypercubic trun- 
cation has only slightly affected the Ginsparg- 
Wilson property. 

6. Parametrization of the QCD FP Dirac 
operator 

To parametrize the QCD FP Dirac operator 
Dpp for massless quarks, we use a general Dirac 
operator as defined in eq.(Q) with all the couplings 
of the 1-hypercube and all elements of the Clifford 
algebra. The gauge invariant functions we use as 
couplings of D-pp are polynomials in local fluctua- 
tions of the gauge fields. Furthermore we include 
the possibility to smear the gauge fields and to 
project them back to SU(3), i.e. we are using so 
called fat links. We use the new parametrization 
of the FP gauge action [|l9j which also makes use 
of fat gauge links. 

The QCD FP action is also quadratic in the 
fermion fields and the Renormalization Group 
step can again be done analytically. The QCD 
FP equation is the generalization of the free equa- 
tion (?]) given by including gauge fields. In case 
Dp has zero modes, the QCD FP equation is most 
conveniently written as 



D C (V) = Kl 

- n 2 n(u)[Dp(u) 



(10) 

«;Ot( J 7)o(C/)]- 1 fit( J 7) ) 



where n is an optimizablc free parameter of the 
block transformation and U and V are the gauge 
fields on the fine and coarse lattice, respectively. 
They are related through the FP equation of the 
pure SU(3) gauge theory 



iFP 



(V) = min (S 



FP 



(U)+T(U,V)), 



(11) 



d + d ] = dt a 



(9) 



where S FP is the FP action of the pure SU(3) 
gauge theory and T(U, V) is the blocking kernel 
of the block transformation. As the fluctuations 
of the gauge fields on the fine lattice are much 
smaller than those of the coarse gauge fields, the 
Dirac operator used on the fine lattice has to have 
good chiral properties on gauge fields with small 
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fluctuations. This however is by far easier than 
to have a Dirac operator with good chiral prop- 
erties on gauge fields with large fluctuations. An 
important fact for the parametrization of the FP 
Dirac operator is that eq.([l(]) can also be given 
in terms of the propagators 



D c 1 (v) = - + n(u)D- 1 (u)nHu), 



(12) 



as long as Dp has no zero modes. In contrast 
to eq.(10) the equation for the propagator gives 
much more weight to the small physical modes of 
the Dirac operator and can therefore be used to 
improve the properties of the small modes of the 
parametrized FP Dirac operator. 

The parametrization is an iterative procedure. 
We first start at a large value of f3, generate 
thermal coarse gauge configurations V with the 
FP gauge action and determine the correspond- 
ing fine configurations U via minimization as in 
eq.([ll]). Using the free FP Dirac operator on 
the fine configurations (which have very small 
fluctuations), Dq(V) and Dq 1 (V) are calculated 
from eqs.(|To|) and (fl2"|). The couplings of the 
parametrized Dirac operator D pa _ r are determined 
by minimizing the following x 2 -function: 



= ^wd^v - d c v\\ 2 + 



X 



\D- 



(13) 



where the sum runs over a number of random 
vectors v on different configurations and A is a 
weighting factor. The use of vectors for the cal- 
culation of a x 2 function for D par is mandatory 
because the definition of Dq requires a matrix 
inversion which we can afford only for a limited 
number of vectors. The minimization of the \ 2 
yields a parametrized FP Dirac operator D paT (V) 
which has good chiral properties over a larger 
range of gauge couplings than the initial trun- 
cated free FP Dirac operator. 

The fitted parametrized operator D pav is now 
used on fine configurations U', determined via 
minimization from coarse configurations V gen- 
erated thermally at a smaller value of (3. Min- 
imizing the x 2 m eq. (|l3|) again gives D par (V) 
which performs well on an even larger range of 
gauge couplings. The whole procedure is repeated 



until we reach f3 = 3.0 which corresponds to 
/^Wilson ~ 5.75. During this procedure we keep 
leading terms in the naive continuum limit fixed 
such that the tree level mass is zero, the 0(a) 
Symanzik condition is fulfilled, the dispersion re- 
lation starts with slope 1 and the normalization 
of the topological charge is correct jL4j. Further- 
more we fix the free field limit such that we re- 
cover the truncated free FP operator on the triv- 
ial gauge configuration. 

In order to parametrize the operator R in eq. (^) 
we proceed in a similar way as for the Dirac opera- 
tor. We also use a general operator with fat links 
and fluctuation polynomials. The parametriza- 
tion of R is however simpler as it is trivial in Dirac 
space and therefore contains a smaller number of 
operators. In contrast to the equation for the 
block transformation of -Dfp the corresponding 
equation for R eq.(j^) contains no inversion and 
therefore the x 2 ~f unc tion which we minimize in 
order to find the optimal parametrization of i? par 
can be defined as 



x 2 = ii J Rp ar -i?cir 



(14) 



where the norm here is the matrix norm \ \A\\^ 

E 



7. Preliminary results 

We have some preliminary results for the 
parametrized QCD FP Dirac operator. In Figure 
^, we show the spectrum of exact eigenvalues of 
d = V2R D V2R at ft = 3.0, i.e. Pwomm ~ 5.75, 
corresponding to a lattice spacing of a « 0.16 fm. 
We show all the eigenvalues on a 4 4 volume 
(crosses) and those closest to the origin on an 8 4 
volume (circles). We see that the spectrum lies 
quite close to the Ginsparg- Wilson circle, with a 
small additive mass renormalization. This is evi- 
dence that the parametrization is sufficiently rich 
to reproduce well the chiral Fixed-Point proper- 
ties. 

In order to measure hadron masses, one needs 
to calculate the quark propagators. The quark 
propagators are determined by solving the equa- 
tion X = D^-q for some source vector r\. We 
have measured the accuracy in calculating the 
inverse, = \\D(D~ 1 r])k — rj\\, after k conju- 
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Figure 4. Exact eigenvalue spectrum of 
the parametrized FP Dirac operator at (3 = 
3.0 (Ayilson ~ 5.75). 



Figure 5. 10 smallest eigenvalues of A* A at 
(3 = 3.4 (/^Wilson ~ 5.95) for the parametrized 
FP (squares) and Wilson (crosses) Dirac opera- 
tors for 30 configurations. 



gate gradient steps, i.e. how many times 77 is 
multiplied by D in calculating D^r). To achieve 
a given accuracy, we find that the parametrized 
FP Dirac operator requires roughly 3 times fewer 
steps than the unpreconditioned Wilson operator 
corresponding to the same effective pion mass for 
the same gauge configuration. This reduces the 
overhead of D par in production runs. 

The chiral properties of a Dirac operator are 
improved by using the Overlap construction. 
Given some input Dirac operator D , a Dirac op- 
erator which satisfies the Ginsparg- Wilson rela- 
tion is given by 

D ov = \--^=, A=l-D . (15) 
VA^A 

If the input operator already satisfies the 
Ginsparg- Wilson relation, then A'A= 1 and 
-Dov = Do- The inverse square root is calcu- 
lated by some expansion in A* A. How quickly the 
expansion converges is determined by the small- 
est and largest eigenvalues of A^A, which tells us 
how far the input operator is from satisfying the 



Ginsparg- Wilson relation. 

In Figure||, we plot the 10 smallest eigenvalues 
of A^A for 30 configurations of size 8 4 at (3 = 3.4, 
i.e. /Jwilson ~ 5.95, using d (squares) and -Dwiison 
(crosses) as input. The largest eigenvalues of A^A 
are ~ 1.5 and ~ 41 for d and Z?wiison respectively. 
The parametrized FP operator d approximately 
satisfies the Ginsparg- Wilson relation and the in- 
verse square root can be determined accurately 
by, for example, a Legendre expansion with only 
a few terms, which would very much improve the 
already good chiral behavior of d. The expansion 
converges faster by projecting out the smallest 
eigenvalues, which are treated exactly. Using d, 
we see A' A has few small eigenvalues and so the 
projection can be done quite cheaply. In contrast, 
using Dwiison to construct A, most of the eigen- 
values of A^ A are far from 1 and very expensive 
expansions with up to hundreds of terms are re- 
quired to determine the inverse square root p0[ . 
There are also more very small eigenvalues to be 
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Figure 6. Histogram of number of multiplications 
to find 10 smallest eigenvalues of A^A for all con- 
figurations plotted in Figure ||. 



projected out to improve the convergence, which 
would also be quite expensive. 

In Figure ^, we plot a histogram of the num- 
ber of A^A multiplications required to find the 10 
smallest eigenvalues of A' A using d and Dwiison 
as input which were plotted in Figure | Using 
d in the Overlap construction, fewer small eigen- 
values need be treated exactly and these can be 
determined with fewer multiplications than are 
necessary using -Dwiison- 

8. Conclusions 

The project of parametrizing the QCD FP 
Dirac operator is ongoing. Preliminary re- 
sults show that the parametrization achieved so 
far gives a Dirac operator which satisfies the 
Ginsparg- Wilson relation quite well with a small 
additive mass renormalization. Although a ma- 
trix vector multiplication with Z? par is roughly 40 
times more expensive than with -Dwiison, we see 
that calculating eigenvalues and propagators con- 
verges much faster. The chiral properties of a 
Dirac operator are greatly improved by using it 



as input for the Overlap construction. As D paI 
is already quite close to satisfying the Ginsparg- 
Wilson relation, the Overlap construction using 
D par is comparable in cost to using -Dwiison- 
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